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SECTION  I 


USER  GENERAL  INSTRUCTIONS 


(1)  The  program  language  for  PC-1500  Pocket  Computer  is  BASIC  and  the 
elementary  part  of  BASIC  program  instructions  are  given  in  the 
SHARP  Instruction  manual,  which  comes  with  PC-1500  Pocket  Computer. 

(2)  The  8  K  RAM  should  be  attached  to  the  computer  to  increase  its 
capacity.  Printer  is  also  required  to  attach  the  computer. 

(3)  The  program  is  called  "STAT".  This  program  name  is  used  to  load  or 
save  the  program  into  cassette  tape  recorder. 

(4)  The  programs  are  started  in  mode  RUN  by  instruction  RUN  (press  the 
keys  [R]  [U]  [N]  and  [ENTER]  ).  User  is  guided  through  the 
programs  by  simple  questions.  The  user  types  the  chosen  answer  and 
presses  the  key  [ENTER]. 

(5)  The  program  is  recorded  on  magnetic  cassette  tape.  To  load  the 
program  from  the  cassette  recorder  to  the  computer,  press  [C]  [L] 

[0]  [A]  [D]  ["]  [S]  [T]  [A]  [T]  ["]  and  [ENTER]. 

(6)  The  STAT-program  contains  mainly  two  sub  programs,  characterized  by 
the  computations  of  (i)  Single-level  data  and  { i i )  Multi-level  data. 

(7)  For  single-level  data,  the  maximum  number  of  data  points  which  can 
be  input  is  48  and  for  multi-level  data,  it  is  the  same  at  each 
level.  The  capacity  for  the  stress  levels  input  is  7. 

(8)  The  input  data  value  has  to  be  positive  number.  Otherwise  the 
computation  is  resulted  in  error,  indicating  that  an  illogical 
calculation  has  been  made. 

(9)  The  Weibull  parameters  are  estimated  by  the  maximum  likelihood 
method. 
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SECTION  II 


t 
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CONTENTS  OF  PROGRAMS 

(A)  Computation  of  Single-level  Data: 

(1)  Data  input 

(2)  Mean  and  standard  deviation 

(3)  Wei  bull  parameter  estimations 

(4)  Confidence  interval  and  A  &  B  allowables 

(5)  Plot  of  Wei  bull  failure  distribution 

(B)  Computation  of  Multi-Level  Data: 

( 1 )  Data  i nput 

(2)  Weibull  parameter  estimation  (pooled) 

(3)  Plot  of  failure  distributions 

(4)  S-N  curve  (Linear  regression) 


SECTION  III 


DESCRIPTION  OF  STATISTICAL  METHODS 


(A)  Single-Level  Data 

(1)  It  is  assumed  that  the  failure  data  X-j,  X 2,...,  X^  of 
sample  size  n  were  generated  from  a  Wei  bull  distribution  with 
parameters  «  and  b  »  that  is, its  cumulative  distribution  function  is 
expressed  as 

* 

F  (x)  =  1  -  exp  -  ( 

» 

This  data  may  be  obtained  from  ultimate  strength  tests  or  fatigue  test  of 
composite  materials  or  metals.  Note  that  a  and  6  are  called  the  shape 
and  scale  parameters  respectively. 

(2)  The  mean  and  standard  deviation  are  computed  simply  by 


mean  x  =  r.x 

v.x2  -  ():x>2/n  1/2 

standard  deviation  s  =  - 

(31  Weibull  parameter  estimates  are  obtained  by  the  maximum 
likelihood  method;  the  scale  parameter  estimate  &  is  the  solution  of 
the  following  M.L.  equation  (See  [3]): 

A 

>:xu  nnx  1  ?:tnx  „ 

— — x -  -  — - —  -  -  ^  U 

>:xa  a  n 
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The  solution  is  calculated  by  using  the  newton's  approximation  method. 
The  scale  parameter  6  is  estimated  by 


(i 


1/. 

a 


(4)  Confidence  interval  for  3  and  A  &  B  allowables  are  obtained 
from  the  following  formulas: 

95%  C.I.  for  6  =8  exp  [-*.95/*  ] 

a 

A  A 

where  a,  3  are  M.L.E.  of  a  and  6 

and  2  95  =  Crit*cal  values  given  in  [3], 

A  -  allowable  *  B*  [sin  ]1/ct 

B  -  allowable  =  B*  [sin  (V^O) 
where  B*  =  95%  C.I.  for  3. 

(5)  Note  that  if  a  random  variable  X  has  a  Wei  bull  distribution 
with  parameter  a  and  8  ,  then  its  log  transformation  T=lnX  has  the 
extreme  value  distribution  with  parameters  a  *  1/  a  and  b  s  In  8  ;  i.e. 
the  cumulative  distribution  of  T  is 

F  (t>-  1-exp  [-  exp  (  ]. 
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The  graphing  of  cumulative  failure  distribution  is  performed  under  (i) 
without  the  log  transformation  or  ( i i )  with  the  log  transformation.  The 
sample  data  points  are  plotted  using  the  median  rank  p .=^"*^/(n+.4). 

J 

(B)  Mult-Level  Data 

(1)  Under  fatigue  stress  level  S. ,  i=l,2,...,m,  the  fatigue  life 

cycles  X.,,  X.„,  ,  X.  are  observed  from  the  life  test,  which  is 

1 »  i  i n-j 

theoratically  generated  from  the  Weibull  distribution  with  parameters  u 
and  3.  .  It  is  also  assumed  that  the  shape  parameter  a  is  the  same 
for  al 1  stress  level s. 

(2)  Weibull  parameters  are  again  estimated  by  the  M.L.  method.  The 
common  shape  parameter  a  is  estimated  by  solving  the  M.L.  equation  (See 
[2]): 


V 

1 


>:x 


ij 


,nxij 


ry  ci 

fio 


L'lnX 

J _ _ 

ni 


=  0 


The  scale  parameter  3^  (under  stress  levels  S. )  is  obtained  from 


6i  = 


1 


Y.  X.°f 


"i  j  ^ 


1/S 


i— 1,2,. ..  ,m 


(3)  The  plotting  of  failure  distribution  curve  is  performed  exactly 
the  same  way  as  illustrated  in  (A)-(5).  Here  only  the  log  transform  of 
the  data  is  used  as  data  points  in  the  plotting. 

(4)  S-N  curve  fitting  (linear  regression)  is  carried  out  under  the 
assumption  of  Weibull  regression  model  (See  [1]).  The  linear  equation  to 
fit  is 


(a)  Y.=ln3..  =  a  +  b  In  S.  (classical  S-N  curve)  or 


(b)  Y .  =  lnBi  =  a  +  b.  . 


The  case  (b)  may  be  applied  to  the  rupture  testing  data. 
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Because  of  various  sample  size  of  different  stress  level,  we  have  to 
use  the  weighted  least  square  method  (See  [4]). 

The  regression  parameters  a  and  b  are  to  estimated  minimizing 

m  (Y.-Y.)2  ... 

Q  =  z  —  y-1  -  where  Y-  =  a  +  b  In  S.  and  W.  =  1/n.  . 
i=l  wi 

The  following  formulas  are  used  in  our  computations: 

m  (Yi-Yi)2 

RSS  =  residual  sum  of  square  =  I  — n — 

i=i  wi 

RMS  =  residual  mean  square  *  RSS/(m-2) 

m  (Yi  -  Y)2 

SSy  =  total  variation  of  y  *  l  - n - 

i=i  wi 

where  Y  =  mean  of  data  Y. 

correlation  coefficient  =  |  1  -  RSS/SSyj  ^ 2 

SS^g  «  Sum  of  square  for  regression  =  SSy  -  P.SS 

F-score  for  regression  *  (m-2)  SS^/RSS. 

which  has  F-distribution  with  (1,  m-2)  degree  of  freedoms. 

For  a  new  stress  level  S,  the  prediction  of  Bis  computed 

/S  /N 

8  (predicted)  *  exp  (a  +  b  In  S)  (Case  (a)). 

or  .  . 

B  (predicted)  =  exp  (a  +  b  S)  (Case  (b)). 
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SECTION  IV 


PROGRAM  DESCRIPTION 


/ A )  Flow  Diagram: 

Single-Level  Data  Multi-Level  Data 
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(B)  Key  Operation  Procedure: 


Following  data  are  used  to  illustrate  the  key  operations 
(i)  Single-Level  Oata:  T300/934  ( 0/45/90/-45 ) s 
Ultimate  Tensile  Strength  in  MPa 


402 

430 

443 

457 

466 

473 

484 

489 

496 

503 

509 

518 

(ii)  Multi-Level  Data:  T300/934  ( 0/45/90/-45 ) $ 


Tension-Tension  Fatigue  at  3  Levels 


1. 

At 

level  1 

of  maximum 

stress 

290  MPa 

1  872 

851 

2 

!  359 

300 

1  437 

800 

1 

698 

100 

215 

984 

864 

420 

2. 

At 

level  2 

!  of  maximum 

stress 

345  MPa 

87 

373 

60 

912 

51 

260 

11 

491 

81 

571 

69 

711 

59 

320 

51 

848 

54 

187 

3. 

At 

level  2 

1  of  maximum 

stress 

400  MPa 

1 

664 

7 

954 

690 

1 

643 

6 

430 

100 

5 

705 

8 

430 

1 

121 

5 

329 

7 

232 

5 

347 

478 
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DISPLAY 


INPUT 


PRINT  OUT  &  REMARKS 


RUN 

[ENTER] 

MULTI-LEVEL  DATA  ?  (Y/N) 

N 

[ENTER 

HOW  MANY  DATA  N  =  ? 

13 

[ENTER] 

SINGE -LLUEL  DATA 

DATA  -  1 

402 

[ENTER] 

INPUT  DATA  N--  ;  3 

407 .  m 

4 .30.  00 

DATA  -  2 

430 

[ENTER] 

443. 00 

477.  00 

466 . 00 

? 

473. 07 

• 

• 

478.00 

484 . 00 

• 

• 

483. 00 

496. 00 

• 

• 

603.00 

509 . 00 

DATA  -  13 

? 

518 

[ENTER] 

5 i 8 . 00 

MEAN  &  ST.  DEV.  ?  (Y/N) 

Y 

[ENTER] 

MEAN--  4  '2  .  .  V 

ST . DEJ. =  93... 

WEIBULL  PARAMETER  ESTIMATE 

ALPHA-  18.787 

BETA--  4 . 8 7 0 7 E  71 

A, 6-ALLOWABLE  ?  (Y/N) 

Y 

[ENTER] 

357  L.  I  .  CGR  BETA-- 
493. • ■ 

A -ALLOWABLE  - 

377.  '->4 

b-allowable  - 

4  ;  :i .  >•'* 

WEIBULL  PLOT  ?  (Y/N) 

Y 

[ENTER] 

f 

LOG  TRANS.  OF  DATA  ?  (Y/N) 

N 

[enter] 

* 

% 

t 

* 

} 

4 

* 

y 

pS-1 
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DISPLAY 

INPUT 

PRINT  OUT  &  REMARKS 

MULTI-LEVEL  DATA  ?  (Y/N) 

Y 

[ENTER] 

NUMBER  OF  LEVELS  = 

3 

[ENTER] 

HUL  T ]  -LEUEL 

DA  i  m 

AT  LEVEL  -  1 

NO. Of  L FUELS  -  8 

NO.  OF  DATA  = 

6 

[ENTER] 

LEUEL  WITH 

LARSES! 

1  r AT  I  HUE  CYfL.  ' 1 

DATA  IN  CYCLES 

LEUEL- 

] 

DATA  -  1 

NO. or  OAT  A 

-  0 

[ENTER] 

187V'  US! 

? 

1872851 

1  43,’BOO 

2!  ‘>884 

DATA  -  2 

2303300 

1638100 

? 

1437800 

[ENTER] 

864420 

STR.  LEUEL- 

230. 0 

• 

• 

NO . OF  DATA 

4 

83338 

• 

1  1481 

63320 

• 

60312 

DATA  -  6 

? 

864420 

[ENTER] 

8138; 

61  848 

61260 

633!  1 

APPLIED  STRESS  LEVEL  = 

290 

[ENTER] 

STR.  LEUEL  • 

64  1  8  7 

34  3. 0 

AT  LEVEL  -  2 

NO. OF  DATA 

12 

1  664 

NO.  OF  DATA  = 

9 

[ENTER] 

1  643 

670:,' 

DATA  -  1 

87373 

[ENTER] 

6323 

?334 

? 

6430 

84  20 

DATA  -  2 

11491 

[ENTER] 

■>3  -13 

630 

? 

!  IW 

1  1  '  1 

334  ?  ' 

• 

• 

• 

STR.  LEUEL- 

400. 0 

• 

DATA  -  9 

54187 

[ENTER] 
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DISPLAY 


INPUT 


APPLIED  STRESS  LEVEL 
AT  LEVEL  -  3 
NO.  OF  DATA  = 

DATA  -  1 

DATA  -  ? 

? 


345 

12 

1664 

1643 


[ENTER] 

[ENTER] 

[ENTER] 

[ENTER] 


DATA  -  12 
? 

APPLIED  STRESS  LEVEL 


5437  [ENTER] 
400  [ENTER] 


DISPLAY 


INPUT 


PRINT  OUT  &  REMARKS 


WEIBULL  PARAMER  ESTIMATE 


WEIBULl  PLOT  ?  (Y/N) 


Y  [ENTER] 


S-N  (BETA)  CURVE  ?  (Y/N) 
LOG  TRANS  OF  S  ?  (Y/N) 


Y  [ENTER] 

Y  [ENTER] 


SELECT  STRESS  LEVEL  = 


250  [ENTER] 


MORE  STRESS  LEVEL  ?  (Y/N) 


N  [ENTER] 


POOLED  ALPHA  .81:, 
BLTA  -1  -  1.H4V3F:  06 
RE 7  A  3-  6 .  ]  6] '  >E  04 
BL  1  A-  '-i  b.  03 


LQGCCYCLE 7 


I 


RE6RE3SJN  LiNl'  j 
^N':  u  +  b  F  I  c1') 
( N-BC7AP 


RE  • .  MEAN  suuare  - 
'  . 0b3AE- J J 

LllK .  ■  i ii  P f  .  ** .  3 j  < 

1  LORE  Euk  k'Ehk. 

1 860. i 

>  ■  I  f  > .  u F-  PRCE.mri- 
c ) ,  )  ) 

b  ‘ SLOPE  i  ■ 

-13.31 

Q  LINE  FRCEL  1 

'  .  I  460E  0. 

STR.  LEUEL=  3b0.0 
Pt  TA  l  'REU  J  P  T  FJ- 
0.0J80E  03 


(C)  Memory  Contents: 


MEMORY 


DESCRIPTION 


MEMORY 


DESCRIPTION 


(Memory  in  General) 


G$ 


L$ 


C(I) 


D(I) 


X(J) 


Y(J) 


Index  variable 


No.  of  levels 


Index  variable 


Index  variable 


Control  variable 


Control  variable 


Control  variable 


Control  variable 


Control  variable 


n.  =  sample  size  of 
ith  level 


S..  =  stress  level 


x  -  coordinate 


y  -  coordinate 


t  expressed  in  term  of  memories 


(Single  level  data) 


Sample  size 


Mean 


Standard  deviation 


95%  C.I.  for  Beta 


B(I) 


(Wei  bull 


A  -  allowable 


B  -  allowable 


Table  values  from  13) 


parameter  estimation; 


Xa  InX 


l  Xa 


I  InX 


I  x“(lnX)‘ 


E,nX/ni) 


*  hear  X  =  failure  data 
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MEMORY 

DESCRIPTION 

MEMORY 

DESCRIPTION 

G  + 

W  +  C/A/A 

L 

Z  n,  6j  S, 

A 

Shape  parameter  a 

M 

Z  n,  6(2 

B(  I ) 

Scale  parameter  8^ 

G 

K*T  -  1*1 

Weibull  plot 

D 

Intercept  a 

N 

Sample  size 

F 

Slope  b 

Z(j) 

Data 

E 

Residual  Sum  of  Square 

w 

Width  of  graph 

V 

Residual  Mean  Square 

E 

Plotting  Interval  in 
x-coordinate 

R 

Correlation  Coefficient 

S 

Scale  Inticator 

S 

F-score  for  regression 

T 

Scale  Inticator 

L 

New  Stress  Level 

(S-N  Cur 

re) 

Y 

Beta  Predicted 

T 

I  n. 

0 

£[ZX“  lr.x/EXa-] 

P 

IB:  B 

W+ 

I  [u/s  -  (R/S)  J 

Q 

I  ni  (Bi-8i)2 

F+ 

0  -  C/A  -  V 

I 

Z  n,  S, 

J 

in.  8. 

K 

__  E  ni  s? 

t  expressed  in  term  of  memories 
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(0)  Sample  Programs: 

1.  T 300/934  [ 0/± 45/90]  ,  Ultimate  Tensile  Strength  in  MPa  and 
Tension-Tension  Fatigue  Test  Data. 


S1NGE-LEUEL  DATA 


INPUT  DATA  N=  13 

402.00 

430. 00 

443.00 

452. 00 

466.00 

473.00 

478.00 

484.00 

489.00 

496.00 

503.00 

509.00 

518.00 

MEAN-'  472.92 
ST . DEU. -  33.22 

ALPHA--  18.782 
BETA---  4.8702E  02 


95V.  C.l.fOR  BETA - 
473.11 

A -ALLOWABLE  - 

370.34 

B-ALLOWABLE  - 

419.69 


multi -leuel  data 

NO. OF  LEUELS  =  3 

LEUEL  WITH  LARGEST 

FATIGUE  CYCLE  IS 
LEUEL- 1 

NO. OF  DATA  =  6 

1872851 

1437800 

215984 

2359300 

1698100 

864420 

STR.  L£UEL=  290.0 

NO.  OF  DATA  •---  9 

873/3 

11491 

59320 

60912 

81571 

51848 

51260 

69711 

54187 

STR.  LEUEL =  345.0 

NO.  OF  DATA  12 

1664 
1643 
5705 
5329 
7954 
6430 
8430 
7 232 
690 
100 
1121 
5347 

STR.  LEUELS  400.0 

POOLED  ALPHAS  1.815 

BETA-1-  J.5473E  06 

BETA-2=  6. 16J5E  04 

BETA -3-  5.0738E  03 


REGRESS  IN  LINE  IS 
i  ogN-  a  +  b  t  1 09S 
<N=BETA) 


RES.  MEAN  SQUARE  ^ 

1 . 0570E-01 

COR. COEFF . -  0.939 

F-SCORE  FOR  REGR.- 
1360. 13 

DEG.  OF  FREEDOM- 

(1,1) 

b  (SLOPE)  5 
-17. 71 

0  (INTERCEPT)  •= 

1  .  1460E  02 

STR.  LEUEL-  250.0 
BETA  PREDICTED- 
2.0180E  07 


i 
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2.  T300/934  [o/±45/9o]  ,  Static  Tension  Strength  in  MPa  and  Tension-  Tension 

Fatigue  Test  Data 


SINGE-LEUEL  DATA 

INPUT  DATA  N -  20 

345.00 
326.00 
501 , 00 
398.00 
399.00 
^94,00 
402.00 
412.00 
494 . 00 
492.00 
412.00 
418.00 
482.00 
428. 00 
425.00 
444.00 
465. 00 
462.00 
444.00 
456.00 

MEAN*  440.20 
S  T , OEU . =  4  4.00 

ALPHA  -  12.3.13 
BETA-  4.5926E  02 


95*  C .  I  FOR  BETA* 
443  82 

A-ALLOuJABLE  = 

305.47 

B-ALLQUABLE  = 

369. 69 


MULT  I -LEUEL  DATA 

NO. OF  LEUELS  =  6 

LEUEL  uJITH  LARGEST 
FATIGUE  CYCLE  IS 
LEUEL- 1 

NO. OF  DATA  -  3 

2826800 

1698100 

864420 


STR. 

LEUEL* 

42.0 

NO.  OF 

DATA 

*  3 

279440 

351900 

564720 

STR. 

LEUEL* 

46.0 

NO.  OF 

DATA 

*  3 

26764 
58530 
51260 

STR. 

LEUEL - 

50.0 

NO.  OF 

DATA 

*  3 

9888 
7798 
30520 

STR. 

LEUEL* 

54.0 

NO.  OF 

DATA 

*  3 

5705 
1643 
1664 

STR. 

LEUEL* 

58.0 

NO.  OF 

DATA 

•*  3 

8 

229 

1079 

STR. 

LEUEL* 

62.0 

POOLED  ALPHA--  !  ,  591 
BETA  -  1-  l , 9226E  06 
BETA -2=  4. 093. IE  05 
BETA- 3=  4.6244E  04 
BETA-4-  1. 7824E  04 
BETA- 5=  3.3384E  03 
BETA-6*  5.6960E  02 


REGRESS  IN  LINE  IS 
I ogN-  a  +  b  *  I09S 
( N*BETA ) 


RES.  MEAN  SQUARE  * 
2. 5099E-0 1 

COR . COEFF . =  0.996 

F  -SCORE  FOR  REGR.* 
534.91 

DEG. OF  FREEDOM* 

(3,4) 

b  (SLOPE)  * 

-20. 54 

0  (INTERCEPT)  ■* 

9. 1403E  01 

STR.  LEUEL-  38.0 

BETA  PREDICTED* 

1 ..  2252E  02 

STR.  LEUEL-  34.0 

BETA  PREDICTED* 

1 . 6956E  08 
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3.  T/300/934,  [02/+452/902]s>  Tension-Tension  Fatigue  Test  Data 


blNOE  LEUEL  DATA 
INPUT  DATA  N-  20 

l 1491,00 

12578.00 
40270.00 
41200. 00 
44830.00 
51848. 00 
54187.00 
58530.00 
59320. 00 
60912.00 
64070.00 
697  t 1  00 
70049.00 
,->04  97.  00 
7. 14 00. 00 
91571  00 
87373  00 
1  1  6667  00 
367644 . 00 
s  i  2600. 00 

DEAN-  97587.40 
9T.OEU.  171686.04 

ALPHA-  1  071 
BETA-  1..0086E  05 


95V.  C.l  FOR  BETA  - 
68085.68 
A -ALLOWABLE  - 

978. 74 

h  allowable  - 

8379. 71 


MULTI -LEUEL  DATA 

NO. OF  LEUELS  =  2 

LEUEL  WITH  LARGEST 

FATIGUE  CYCLE  IS 
LEUEL -1 

NO.  OF  DATA  20 

1  1491 

17578 

40270 

4  1200 

4  4  830 

51848 

54187 

58530 

59320 

60912 

64070 

69711 

70049 

7049  7 

7  1400 

8  1571 
6  7373 

1  .1  6667 
367644 
5  >  2600 

STR.  LEUEL-  50.0 

NO.  OF  DATA  •-  20 

61 
100 
359 
578 
690 
1121 
1588 
i  64  T 
1  064 
<558 
4  930 
329 

1  7 
6  j5 
6430 
7232 
7954 
8430 
9060 
9590 

STR.  LEUEl-  58.0 


POOLED  ALPHA--' l  .  035 
BETA  - 1  -  9.9J95E  04 
BETA-2-  4. 1 172E  03 


REGRESS  IN  LINE  IS 
1 04 N -  a  +  b  *  I 09S 
(N-BETA) 


b  (SLOPE) 

-2. 1437E  01 
j  (INTERCEPT;  - 
9.5369E  01 

STR.  LEUEL=  42.0 

BF  TA  PRF  0 1 C TED  - 
4  .  1 664 E  06 

STR.  LEUEL-  34.0 

BF  TA  PRLDirTLD- 
'•) .  864  71  08 
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(E)  Program  Listing: 


l: "STAT". TEXT  . 
COLOR  0 

10: DIM  XC50), YC50 
),  ZC50),  AC7,  50 
>,  B<50),  CC7 ),  D 
(?) 

.15:  INPUT  "NULTI-L 
EUEL  DATA  ?  (V 
-'N) "  i  G$ 

20: IF  G$=“V"GOTO 
600 

25: LPR I  NT  "SINGE- 
LFUEL  DATA" . LF 
(  1 ) 

28 • GOSUB  1190 

30: INPUT  "HOW  DAN 
V  DATA  N=?";CC 
1  ) 

35: LPR I  NT  "INPUT 
DATA  N="iC(l): 
LF  (1) 

45:  C  =  l,  K- 1 ,  Z=\ 

50: FOR  J=1T0  CC1) 

60: PAUSE  "DATA  "  + 
STR$  CJ)+"  •=" : 
INPUT  A(l,  J) 

65:  GOSUB  1200 

70:  LPRINT  AC  1,  J): 
NEXT  J 

150:  INPUT  "DEAN  & 
ST.DEU.  ?  (Y/N 
> " ; G$ : IF  G$-"N 
"GOTO  700 
155-’  N-CC 1  ) :  LF  Cl) 

1 70:  T -0,  S=0:  FOR  J  - 
1  TO  N 

174. T -T  fA C 1,  J),  S-S 
t-AC  1,  J  )*A(  1,  J  .) 
136: NEXT  J 
1 38-  S=S~T*T/N 
180:  T=T/N,  S=S/CN-1 
) 

182:  S=TS 
183: GOSUB  1200 
184: LPRINT  "MEAN=" 
;  t 

186: LPRINT  "ST.DEU 
. =" } S: LF  Cl) 
200: GOTO  700 
321 :  IF  K>  1  GOTO  .351 
325- GRAPH  -COLOR  1 
330. GLCURSOR  C25, - 
240). SORGN 


335- LINE  C-5,  0)-Cl 
90,  0),  0,  1 

337: LINE  C 130,  0)-C 
190,  200),  0,  1 
338: L  3NE  C190, 200) 
-C0,  200),  0,  1 
340: LINE  C0, 200)- C 
0,  -5),  0,  1 
345: W=1 80/ZCN) 

346: LINE  C90,  0)-C9 

0,  — 5),  0,  1 

347: LINE  C180, 0)-C 
180,  -5),  0,  1 
348: L 1 NE  C0, 100)-C 
-5,  100),  0,  1 
349; LINE  C0, 200)-C 
-5,  200),  0,  1 
351 :FOR  J  - 1 T 0  N 
352 : LET  YCJ)=CJ-.3 
VCN+.  4 ) ;  NEXT 
J 

360: FOR  J - 1 T 0  N 
362: XCJ)-ZCJ)*W 
364: VCJ)=YCJ)*200; 
NEXT  J 

370: FOR  J-1TO  N 
375: LINE  CXCJ)-2, V 
(J)-2)-(XCJ)+2 
,  YCJ)42),  0,  1 
376: LINE  CXC J )-2,  v 
<J)+2)-(XCJ>+2 
,  VCJ)-2),  0,  l. 
NEXT  J 

385; IF  K > ] GOTO  400 
330: E=ZCN)/40 
400: IF  D$="N"GOTO 
480 

415: FOR  J-0TO  42 
420: XC J )-E*J 
430:  VC  J  )d  -EXP  C- 
EXP  C  CXC J )-LN 
BCJ))*A)) 

440:  XCJ)-X(J) *14,  YC 
J)=YCJ)*200 
450: NEXT  J:GOTO  49 
0 

480: FOR  J  -0TO  42 
485: XC J )-E*J 
486: VC J )- 1-EXP  C-C 
XCJ )/B( I )  )-^A) 
487: XC  J)-XC  J)*LI, VC 
J)-YCJ)*200: 
NEXT  J 


490. FOR  J  —  1 T 0  42 
492: LINE  CXC J -  1 ), V 
C  J~ 1 ) )- CXC J ),  v 
(J)),  0,  1 :  NEXT 
J 

496: IF  Z-2AND  K=1 
GOTO  555 

497: IF  K^CGOTO  500 
498:  IF  K  >  1  GOTO  555 
500: GLCURSOR  C-25, 
195) 

510: LPRINT  "1." 
520: GLCURSOR  C-25, 
95) 

530: LPRINT  ".5" 
540: GLCURSOR  C 50,  - 
20) 

545: S=INT  CEX400) 
550: S=S/10 
551 : T=S/2 

552 : LPR  J  NT  STRS  CT 

) 

553: GLCURSOR  C 1 40, 
-20) 

554 : LPRINT  STR$  (S 
) 

555: RETURN 
558: LF  C2) 

560: INPUT  "A, B-ALL 
0 LIABLE  ?  (Y/N) 
"  j  G$ 

562: IF  G$-"N"GOIO 
900 

564 .GOSUB  1250 
565: GOSUB  1200 
566: H-BC 1 )¥EXP  (-B 
(N)/A) 

567:  LPRINT  "95>:  C. 

I.  FOR  BETA-",  H 
568:  P~H*C CLN  Cl/. 9 
9) )AC 1/A) ) 

570: Q  =  H*C  CLN  Cl/. 9 
) )AC 1/A) ) 

575: LPRINT  "A-ALLO 
LIABLE  -  ",  P 
580: LPR 1  NT  "B-ALLC 
LIABLE  =  ",  Q 
590: GOTO  900 
600: Z-2 : LF  C 2) 

601 : LPRINT  "NULTI  - 
LEUEL  DATA" 
602: LF  Cl) 

603: INPUT  "NUMBER 
OF  LEUELS  =" ; C 
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605- GOSUB  1190 
606: LPRINT  "NO. OF 
LEUELS  ='* » C : LF 
Cl) 

608- LPRINT  "LEVEL 
WITH  LARGEST  F 
ATIGUE  CYCLE  1 
S  LEUEL-1 

;  LF  (1) 

610:FOR  I  =  1  TO  C 
615: PAUSE  "AT  L EVE 
L-"+STR$  (1) 

620: INPUT  "NO. OF  D 
ATA  =  ";C(D 
624: GOSUB  1190 
625 : LPR 1NT  "NO. OF 
DATA  =HjC(]> 

622: (F  I >1G0T0  630 
628: PAUSE  "  DATA 
IN  CYCLES" 

630:  FOR  J-ITO  CCD 
640: PAUSE  "DATA - "+ 
STR$  (J): INPUT 
ACI,  J  ) 

645: GOSUB  1210 
650: LPR1NT  ACI,  J): 
NEXT  J 

663: GOSUB  1192 
665: INPUT  "APPLIED 
STRESS  IEUEL= 

"  ;  DC  1  ) 

625: LPRINT  "STR.  L 
EUE‘L  =  "  i  DC  D  •'  LF 
<  1  >• 

690: NEXT  J 
200  :A---2:  BEEP  2 
705: PAUSE  "HE  1  BULL 
PARANETER  EST 
I  NATES' 

710.U=0.  W=0,  0-0 
720:  FOR  I--ITO  C 
725- R=0,  S=0,  T=0,  U= 
0 

730- FOR  J'-ITO  CCD 
740. R=R+CAC 1,  J  >rtA) 
*LN  ACI, J> 

750. S=S+A( 1,  J>*A 
760- TsT+LN  AC  1,  J) 
770. 0=U+ (AC 1 ,  J>AA> 
*LN  ACI,  J>*LN 
ACI,  J> 

775. NEXT  J 
805: W=W+(U/S-  CR/S) 
*2)> U=U+DC( 1 ) 
8 1 0 . 0=0  *■  C R'S >  •  NEXT 
1 

825.  r--0  -C^A  -U,  G=1«H 
C/CA*A> 


830-  IF  ABS  (F/GK. 

01 GOTO  847 
845: A=A-F/G: GOTO  7 
10 

847: GOSUB  1185 
849;  IF  2 - 1G0T0  852 
850: LPRINT  "POOLED 

alpha=".lf  c- 
1) 

851 : LPRINT  A: GOTO 

£53 

852:  LPRINT  " ALPHA - 
";A 

853; FOR  1*1  TO  c 

855: BCI )=0:NEXT  I 
860: FOR  1  =  1  TO  C 
863:  FOR  J=1T0  CCD 
865:BCD=BCD  +  CAC1 
, J)^A):NEXT  J 
866:BC1)=CBC1)/CC1 
1  )<^C  17A) 

868: GOSUB  1170 
870 ;  IF  2--1  GOTO  872 
871 .LPRINT  "BETA-" 
4-STR*  (1)+"=", 
BCD  .’GOTO  875 
872: LPRINT  "BETA=" 
;BCD 

874:  IF  2=1  GOTO  558 
875: NEXT  1 ;  LF  CD 
900: INPUT  "WEI BULL 
PLOT  ?  CY/N)" 
;G$: IF  G*="N" 
GOTO  1300 
901 : IF  2=2LET  D*=" 
V" : GOTO  904 
902: INPUT  "LOG  TRA 
NS. OF  DATA  ?  C 
Y/N)"iD$ 

904: FOR  J -1  TO  C 
905:N=CCD,  K=J 
906:  FOR  J--1T0  N 
907; 2CJ)=AC  I, J): 
NEXT  J 

909: IF  D$="V"GOSUB 
1260 

910:  GOSUB  1005 
928.  GOSUB  321 
940:  NEXT  1 
94 5 .COLOR  0 • TEXT 
950:  LF  CD 
955: IF  D$="N"GOTQ 
965 

960. IF  2=2G0TQ  962 
961: LPRINT  "  L 

OGCDATA)". GOTO 
965 

962: LPRINT  ”  LO 

GC CYCLE)" 

965: IF  2=1G0T0  15 


970: LF  ( 2 )  •  GOTO  13 
00 

1000: Z(  J  >=Z(H), 2C 
H)=N. RETURN 
1005:  FOR  J--2T0  N 
U10.AREAD  Z  ( J  > 
1120: FOR  H -  1  TO  J- 
1 

1 130.  IF  2CJK2CH) 
LET  f1=Z(J>. 
GOSUB  1000 
11 35: NEXT  H 
11 60: NEXT  J: 

RETURN 

1 170 : US 1NG  "4.444 
4*" : RETURN 
l 185: USING  "444.4 
4# " . RETURN 
l 190: USING  '444". 
RETURN 

1197- USING  "44444 
4": RE  TURN 
1200-USING  "44444 
44.44": 
RETURN 

121 0 : US J NG  "44444 
444". RETURN 
1250: BC5>  =  1 .  102,  B 
(6>  =  .939,  BC7 

>  =  .829,  BC8>  = 

. 751,  BC9)  =  . 6 
91,  B(10>=.64 
4,  BCD  >  =  .605 

1251 : BC 12)=. 572, B 
C 13)=. 544, BC 
14)=. 520, BC 1 
5)=. 499, BC 16 

>  =  . 48, BC 17)= 
.463,  BC  18)  =  . 
447 

1252: BC 19)  =  . 433,  B 
(20)  =  . 421.  BC 
21)=. 409, B ( 2 
2)= . 398, BC23 

>  =  . 388, B ( 2 4  > 
=  .379,  BC25)  = 

37 

1253:8(26)  =  . 362,  B 
(27)  =  . 354,  BC 
28)=. 347, BC2 
9)  =  . 34,  BC30) 
=  .  334,  B ( 3  1  )  = 
328, BC  32 )  =  . 
323 


1254: B<33)=.  317,  B 
<34)-. 312,  B( 
35)-. 307, B(3 
6)-. 302,  B<37 
>  =  . 297,  B<38) 
=  .  293 

1255: B<39)-. 289,  B 
<40)-. 285,  B< 
41)=. 281, B<4 
2)  =  .278,  B<43 
.  274,  B <44 ) 
=  .  271 

1256: B<45)-. 267,  B 
<46)-. 264,  B< 
47)  =  . 261,  B<4 
8)-. 258, B<49 
)  =  . 255,  B  <50) 
=. 253: RETURN 
1260: FOR  J-l TO  N 
1262: Z< J )-LN  Z<J) 

: NEXT  J 
1265: RETURN 
1300: IF  Z=1 GOTO  1 
5 

1302: INPUT  "S-N(B 
ETA)  CURUE  ? 
<Y/N)";L$ 
1303: IF  L$="N" 
GOTO  15 

1305: INPUT  "LOG  T 
RANS.OF  S  ? 
<Y/N)" j  S* 
1306:  IF  S$="N" 
GOTO  1308 
1307: LPRINT  "REGR 
ESSIN  LINE  I 
S  logN-  a  + 
b  ¥  I09S 
<N=BETA) " : 
GOTO  1310 
1308: LPRINT  "REGR 
ESSIN  LINE  I 
S  logN  -  a  + 
b  *  S 
<N=BETA) " 
1310: T-0,  P-0; LF  < 

1  ) 

1312: FOR  B-1TO  C 
1315: IF  S$="N" 
GOTO  1320 
1317: X<B)=LN  D<B) 

:  GOTO  1324 
1320: X<B)=D<B) 
1324: Z<B)=LN  <B<B 
)) 

1326: T-T+C<B):P=P 
+Z(B) : NEXT  B 


1328: 1=0: J-0: K-0: 

L-0: 11=0: Q-0 
1330: P-P/C 
1332: FOR  B-1TO  C 
1335: Q=Q+<Z<B)-P) 
*<Z<B)-P)*C< 
B) 

1340: I=I+C(B)*X(B 
):J=J+Z<B)*C 
<B) 

1350:  K=K+X<B)¥X(B 
)*C<B):L=L+X 
<B)*Z<B)*C<B 
) 

1360:  f1=t1+Z<B)¥Z<B 
)*C<B) : NEXT 
B: LF  <1) 
1410:  G=tOkT-I*I 
1420:  0- <J*K-I*L)/ 
G 

1430:F=<T*L-I*J)/ 

G 

1440: E-0 

1450: FOR  B-l TO  C 
1460: Y(B)=Z<B)-<D 
+F*X<B)) 

1470: Y<B)=Y<B)*Y< 
B)*C<B) 

1480: E-E+Y<B) : 
NEXT  B 

1495: IF  C-2G0T0  1 
540 

1500: U =E/<C-2) 
1505: GOSUB  1170 
1510: LPRINT  "RES. 
MEAN  SQUARE 
",'U:LF 

<1) 

1520: GOSUB  1185 
1522: R-T  < l-E^Q) 
1523: LPRINT  "COR. 
COEFF . =" ; R: 
LF  <  1 ) 

1525: GOSUB  1200 
1526: S=<Q-E)xU 
1527: LPRINT  "F-SC 
ORE  FOR  REGR 
=  11  •  s 

1530: LPRINT  "DEG. 
OF  FREEDOM- 
<1. 

"+STR$  <C-2) 
+")":LF  <1) 
1540: LPRINT  "  b  < 
SLOPE)  = 

"JF 

1543: GOSUB  1170 


1545: LPRINT  "  a  < 
INTERCEPT)  - 
"ID: 

LF  <1) 

1560: GOSUB  1197 
1570: INPUT  "SELET 
STRESS  LEUE 
L  =  "  ,*  0 

1575: LPRINT  "STR. 

LEUEL-" ; 0 
1577: IF  S*="N" 
GOTO  1580 
1578: O-LN  0 
1580: Z=D+F*0 
1 590: Y-EXP  <Z) 
lu95: GOSUB  1170 
1600: LPRINT  "BETA 
PREDICTED- 

M  •  y 

1610: INPUT* "MORE 
STRESS  LEUEL 
?  <Y/N)";R$ 

: LF  <1) 

1620: IF  R$="Y" 
GOTO  1560 
1650: GOTO  15 
1700: END 
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SECTION  V 


CONCLUSION 

The  description  and  instruction  of  the  use  of  Sharp  PC-1500  Pocket 
Computer  for  the  key  calculation  of  statistical  data  analysis,  in  Weibull 
distribution,  are  presented  in  this  paper.  Instant  calculation  of  A  and 
B  allowables,  S-N  curve  fitting  and  prediction  can  be  made  for  practical 
use. 
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